Dynamical localization simulated on a few qubits quantum computer 
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We show that a quantum computer operating with a small number of qubits can simulate the 
dynamical localization of classical chaos in a system described by the quantum sawtooth map model. 
The dynamics of the system is computed efficiently up to a time t > I, and then the localization 
length £ can be obtained with accuracy v by means of order 1/v 2 computer runs, followed by coarse 
grained projective measurements on the computational basis. We also show that in the presence 
of static imperfections a reliable computation of the localization length is possible without error 
correction up to an imperfection threshold which drops polynomially with the number of qubits. 

PACS numbers: 03.67.Lx, 05.45.Mt, 24.10.Cn 



I. INTRODUCTION 



Recent experimental progress in nuclear magnetic res- 
onance (NMR)-based quantum processors allowed the 
demonstration of quantum algorithms [1], including 
Grover's algorithm [2] and quantum Fourier transform 
[3]. More recently, it has been possible to implement the 
simplest instance of Shor's algorithm, namely the fac- 
torization of 15, using 7 qubits and a sequence of about 
300 spin-selective radio- frequency pulses [4]. In parallel, 
thanks to the development of techniques for the manip- 
ulation of cold atoms in linear traps, the realization of 
up to 50 two-qubit control-not gates within the relevant 
decoherence time scale is currently becoming possible [5] . 
Solid state realizations are also under way in several ex- 
perimental groups working with various solid state de- 
vices. In particular, it has been demonstrated that a 
superconducting tunnel junction circuit can behave as 
an artificial spin 1/2 atom, whose evolution can be con- 
trolled by applying microwave pulses. The quality factor 
of quantum coherence is sufficiently high to envisage the 
realization of two-qubit gates based on capacitively cou- 
pled circuits of this type [6] . 

In this context, it is of primary importance to find 
efficient quantum algorithms, that could be usefully sim- 
ulated with a small number of qubits. Such algorithms 
would naturally become the ideal software for demonstra- 
tive experiments in the coming generation of quantum 
processors. Dynamical models represent a natural testing 
ground for quantum information processors. The algo- 
rithm for the quantum baker's map [7] has been recently 
implemented on a three qubit NMR quantum processor 
[8]. These experiments tested the sensitivity to pertur- 
bations, in a system which is characterized by chaotic 
unpredictable dynamics in the classical limit. 

In this paper we show that quantum computers can 
simulate efficiently the quantum localization of classi- 



cal chaos. Dynamical localization is one of the most 
interesting phenomena that characterize the quantum 
behavior of classically chaotic systems: Quantum inter- 
ference effects suppress chaotic diffusion in momentum, 
leading to exponentially localized wave functions. This 
phenomenon was first found and studied in the quan- 
tum kicked rotator model [9] and has profound analogies 
with Anderson localization of electronic transport in dis- 
ordered materials [10]. Dynamical localization has been 
observed experimentally in the microwave ionization of 
Rydberg atoms [11] and is now actively studied in exper- 
iments with cold atoms [12]. 

In this paper, we study dynamical localization for the 
quantum sawtooth map, using the algorithm developed 
in Ref. [13]. This algorithm has some specific advantages 
with respect to similar algorithms for the simulation of 
other dynamical systems, for instance the kicked rotator 
[14]. There are no extra work space qubits, namely all the 
qubits are used to simulate the dynamics of the system. 
This implies that less than 40 qubits would be sufficient 
to make simulations inaccessible to present day super- 
computers. We note that this figure has to be compared 
with the more than 1000 qubits required to the Shor's 
algorithm to outperform classical computations. We will 
also show that in this model dynamical localization could 
be observed already with 6 qubits. 

The paper is organized as follows. In Section II we de- 
scribe the sawtooth map model, focusing on the regime 
of dynamical localization. In Section III we show that 
a quantum computer operating with few qubits can in- 
deed perform simulations of dynamical localization. In 
Section IV we discuss how to extract information (the 
localization length) from the quantum computer wave 
function. In Section V we study the stability of those 
computations in the presence of static imperfections in 
the quantum computer hardware. In Section VI we dis- 
cuss the transition to quantum chaos, induced by static 
imperfections, in the quasienergy spectral statistics. In 
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Section VII we present our conclusions. 



t* « D n w (7r 2 /3)fc 2 , 



(3) 



II. THE MODEL 

The quantum sawtooth map is the quantized version 
of the classical sawtooth map, which is given by 

n = n + k(8 -tt), 8 = 8 + Tn, (1) 

where (n, 0) are conjugated action-angle variables (0 < 
8 < 27r), and the over bars denote the variables af- 
ter one map iteration. Introducing the rescaled mo- 
mentum variable p = Tn, one can see that the clas- 
sical dynamics depends only on the single parameter 
K = kT. The map (1) can be studied on the cylinder 
(p G (— oo, +oo)), which can also be closed to form a torus 
of length 2irL, where L is an integer. For K > 0, the mo- 
tion is completely chaotic and exhibits normal diffusion: 
< (Ap) 2 >ss D(K)t, where t is the discrete time mea- 
sured in units of map iterations and the average < • • • > 
is performed over an ensemble of particles with initial 
momentum po and random phases < 8 < 2tt. For 
K > 1 the diffusion coefficient is well approximated by 
the random phase approximation, D(K) w (n 2 /3)K 2 . 

The quantum evolution in one map iteration is de- 
scribed by a unitary operator U (called Floquet operator) 
acting on the wave function if>: 

$ = Uil> = e- 4Tfi2/ V fe( <^ )2/2 V, (2) 

where n = —id/ 89 and ip(8 + 2tt) = ip(8) (we set h = 1). 
The classical limit corresponds to k — > oo, T — > 0, and 
K = kT = const. In Refs. [13,15,16] we studied the 
map (2) in the semiclassical regime. This is possible by 
increasing the number of qubits n q = log 2 N (N is the 
total number of levels), with T = 2wL/N, K — const. 
In this way, the number of levels inside the "unit cell" 
— n < p < t: {L — 1) grows exponentially with the 
number of qubits (-N/2 < n < N/2), and the effective 
Planck constant h c g ~ h/k <~ 1/N — > when N —* oo. 

Differently from previous studies, in this paper we 
study the map (2) in the deep quantum regime of dynam- 
ical localization. For this purpose, we keep k, K constant. 
Thus the effective Planck constant is fixed and the num- 
ber of cells L grows exponentially with the number of 
qubits (L = TN/2n). In this case, one studies the quan- 
tum sawtooth map on the cylinder (n e (— oo,+oo)), 
which is cut-off to a finite number of cells due to the fi- 
nite quantum (or classical) computer memory. We stress 
again that, since in a quantum computer the memory ca- 
pabilities grow exponentially with the number of qubits, 
already with less than 40 qubits one could make simula- 
tions of systems inaccessible for today's supercomputers. 
Similar to other models of quantum chaos [9], the quan- 
tum interference in the sawtooth map leads to suppres- 
sion of classical chaotic diffusion after a break time 



where D n is the classical diffusion coefficient, measured 
in number of levels (< (An) 2 D n t). For t > t* 
only An <~ D n levels are populated and the localization 
length I ~ An for the average probability distribution is 
approximately equal [17]: 

I ~ D n (4) 

Thus the quantum localization can be detected if I is 
smaller than the system size N. 

In the following we consider K = y/2, two values of 
k, k = \/3 and k = 2, and 6 < n q < 21, so that the 
above analytical estimate gives £(k = \/3) w 10 and 
£(k = 2) w 13 < N. We assume that at t = the system 
is in a momentum eigenstate, ip(n) = 5 nno . Since this is a 
quantum register state, it can be obtained in 0(n q ) onc- 
qubit operations starting from the fiducial state ( "ground 
state") n of the quantum computer hardware. 



III. SIMULATION OF DYNAMICAL 
LOCALIZATION 

An exponentially efficient quantum algorithm for the 
simulation of the map (2) was described in Refs. [13,16]. 
It is based on the forward/backward quantum Fourier 
transform [1] between momentum and coordinate bases: 
On the whole, it requires n g — 3n 2 + n q gates per 
map iteration (3n 2 — n q controlled-phase shifts and 2n q 
Hadamard gates). This number has to be compared with 
the 0(N log N) operations required to a classical com- 
puter to simulate one map iteration by means of the fast 
Fourier transform. 

In Fig. 1, we show that, using our quantum algorithm, 
exponential localization can be clearly seen already with 
n q = 6 qubits. After the break time t*, the probability 
distribution over the momentum eigenbasis decays expo- 
nentially, 

W n = mn)\^\^(-^^), (5) 

with no = the initial momentum value. Here the lo- 
calization length is Ik, 12, and classical diffusion is sup- 
pressed after a break time t* w £, in agreement with the 
estimates (3)- (4). This requires a number N g w 3n 2 £ ~ 
10 3 of one- or two-qubit quantum gates. The full curve of 
Fig. 1 shows that an exponentially localized distribution 
indeed appears at t w t*. Such a distribution is frozen 
in time, apart from quantum fluctuations, which we par- 
tially smooth out by averaging over a few map steps. The 
localization can be seen by the comparison of the proba- 
bility distributions taken immediately after t* (full curve 
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in Fig. 1) and at a much larger time t — 300 w 25t* 
(dashed curve in the same figure). 

We also note that the asymptotic tails of the wave 
functions decay as a power law(see Fig. 2), 



W n oc 



1 



n - no 



(6) 



This happens due to the discontinuity in the kicking 
force of Eq. (1), f(0) = k(6 — n), when the angle vari- 
able 8 — 0. For that reason the matrix elements of the 
one period evolution operator U for the quantum map 
(2) decay as a power law in the momentum eigenbasis: 
U nm = (n\U\m) ~ l/|n — m| Q , with a = 2. This case 
was investigated for random matrices, where it was shown 
that eigenfunctions are also algebraically localized with 
the same exponent a [18]. We also note that dynamical 
localization in discontinuous systems was studied in [19]. 
Since the localization picture is not very sensitive to the 
behavior of the tails of the wave function, a rough esti- 
mate of the crossover between the exponential decay (5) 
and the power law decay (6) is given by their crossing 
point, 



W n (n c ) 



1 



dog< 



(7) 



This implies that by increasing I the exponential local- 
ization is pushed to larger momentum windows and down 
to smaller probabilities. 




FIG. 1. Exact quantum computation of probability dis- 
tribution over the momentum basis with n q = 6 qubits for 
k = \/3 and initial momentum no = 0; the average is taken 
in the intervals 10 < t < 20 (full curve) and 290 < t < 300 
(dashed curve). The straight line fit, W„ oc exp(—2\n\/£), 
gives a localization length I « 12. Here and in the following 
figures, the logarithms are decimal. 



IV. MEASUREMENTS 

We now discuss how it would be possible to extract 
information (the value of the localization length) from 
a quantum computer simulating the above described dy- 
namics. The localization length can be measured by run- 
ning the algorithm several times up to a time t > t*. 
Each run is followed by a standard projective measure- 
ment on the computational (momentum) basis. The out- 
comes of the measurements can be stored in histogram 
bins of width Sn oc £, and then the localization length 
can be extracted from a fit of the exponential decay of 
this coarse-grained distribution over the momentum ba- 
sis. In this way the localization length can be obtained 
with accuracy v after the order of 1/v 2 computer runs. 
It is important to note that it is sufficient to perform a 
coarse grained measurement to generate a coarse grained 
distribution. This means that it will be sufficient to mea- 
sure the most significant qubits, and ignore those that 
would give a measurement accuracy below the coarse 
graining Sn. Thus, the number of runs and measure- 
ments is independent of I. However, it is necessary to 
make about t* ~ I map iterations to obtain the local- 
ized distribution (see Eqs. (3,4)). This is true both for 
the present quantum algorithm and for classical com- 
putation. This implies that a classical computer needs 
0(^ 2 log£) operations to extract the localization length, 
while a quantum computer would require 0(^(log£) 2 ) el- 
ementary gates (classically one can use a basis size N ~ I 
to detect localization). In this sense, for I ~ TV = 2 rlq the 
quantum computer gives a square root speed up if both 
classical and quantum computers perform O(N) map it- 
erations. However, for a fixed number of iterations t the 
quantum computation gives an exponential gain. For 
I <C N such a gain can be very important for more com- 
plex physical models, in order to check if the system is 
truly localized [20]. 



V. EFFECTS OF STATIC IMPERFECTIONS 

In order to study the effects of static imperfections 
on the stability of the above described algorithm, we 
model the quantum computer hardware as a linear ar- 
ray of qubits with static imperfections, represented by 
fluctuations in the individual qubit energies and resid- 
ual short-range inter-qubit couplings [21]. The model is 
described by the following many-body Hamiltonian: 



H s = £(A 



(8) 



i<3 



where the crfi's are the Pauli matrices for the qubit 
i, and A is the average level spacing for one qubit. 
The second sum in (8) runs over nearest-neighbor qubit 
pairs, zero boundary conditions are applied, and 5i, Jij 
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are randomly and uniformly distributed in the intervals 
[—3/2,5/2] and [—J, J], respectively. We model the im- 
plementation of the above algorithm on this hardware 
architecture as a sequence of instantaneous and perfect 
one- and two-qubit gates, separated by a time interval 
T g . Therefore we study numerically the evolution in time 
of the quantum computer wave function in the presence 
of the following many-body Hamiltonian: 

H(r) = H s + H g (r), (9) 

where 

H g (T)=Y,5(T-kT g )h k . (10) 

k 

Here hk realizes the fc-th elementary gate according to 
the sequence prescribed by the algorithm. We assume 
that the phase accumulation given by Ao is eliminated 
by standard spin echo techniques [1]. In this case, the re- 
maining terms in the static Hamiltonian (8) can be seen 
as residual terms after imperfect spin echoes and give 
unwanted phase rotations and qubit couplings. 

The effect of static imperfections on the probability 
distribution over the momentum basis is shown in Fig. 2, 
for k = 2, n q = 11, t = 100, J = 0, and different rescaled 
imperfection strengths e = 5r g . For e = 10~ 4 , the local- 
ization peak is reproduced with high fidelity, while the 
tails of the wave function are strongly enhanced. This 
is due to the fact that errors affecting the most signif- 
icant qubits can induce a direct transfer of probability 
very far in the momentum basis [22,13]. For e = 10~ 3 a 
measurement of the decay of the localization peak would 
overestimate the localization length by a factor two, while 
for e = 10~ 2 any trace of dynamical localization has been 
destroyed. 




-1000 -500 500 1000 



n 

FIG. 2. Probability distributions for k — 2, n q = 11, 
no = 0, J = 0. From bottom to top: e = 0, 1CP 4 (gray 
line), 10~ 3 (shifted up by a factor 4), and 10~ 2 (shifted up 
by a factor 8). The straight line gives a localization length 
15. 



In order to study in a more quantitative way the sta- 
bility of quantum computation in the presence of static 
imperfections, we consider the following two quantities: 

(i) the diffusion coefficient D n (e), obtained from the re- 
lation 

< (An) 2 D n (e)t. (11) 

This is an important characteristic related to transport 
properties of the system; 

(ii) the inverse participation ratio 



this quantity determines the number of basis states sig- 
nificantly populated by the wave function and gives an 
estimate of the localization length of the system. We 
stress that, differently from the previous quantity, £ is 
local in the localized regime, i.e. it is insensitive to the 
behavior of exponentially small tails. 




5 10 15 20 



FIG. 3. Dependence of the wave function second moment 
on time, for k = 2, n q = 11, J = (full symbols) and J = 8 
(empty symbols) at e = (circles), 5 x 10~ 5 (squares), 10~ 4 
(diamonds), and 2 x 10 -4 (triangles). The straight lines fits 
give the diffusion coefficient D n (e). The curves are aver- 
aged over 10 disorder realizations and 10 initial conditions 
n e [—5, 5[. 

In Fig. 3 we show < (An) 2 > as a function of time, 
for n q = 11 qubits, J = 0, and different imperfection 
strengths e. By means of these curves we extract the 
diffusion coefficients D n (e) from linear fits extended to 
the first few map steps. In the same figure, we show 
that similar curves are obtained for J = S. The de- 
pendence of the inverse participation ratio on t is shown 
in Fig. 4, again for /c = 2, n q = 11. We note that, 
for imperfection strengths strong enough to induce huge 
variations in the diffusion coefficient (D n (e) 3> D n (0j), 
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£ is only slightly modified (£(e) ~ £(0)). Iterating the 
map (2) long enough (t > £(e)), we get the saturation 
value £00 (e)- This quantity increases with e and one has 
complete derealization when ^ (e) ~ N (this is evident 
for e = 5 x 10~3 in Fig. 4; in this case £ saturates af- 
ter t < 100 map iterations). Again we note that similar 
curves are obtained for J = 6 (see Fig. 4). 



in good agreement with the data of Fig. 7. 
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FIG. 4. Time dependence of the inverse participation ratio 
£ for k — 2, n q = 11 qubits, J — (solid lines) and J = 5 
(dashed lines). From bottom to top: e = 0, 5 x 1CT 4 , 1CF 3 , 
2 x 10~ 3 , 5 x 1CT 3 . The curves are averaged as in Fig. 3. 

In Fig. 5 we plot the dependence of diffusion co- 
efficient D n on e for different n q values. From each 
curve we extract the critical imperfection strength er> (n q ) 
corresponding to doubling of the diffusion coefficient, 
D n (e D ) = 2D n (0). The data of Fig. 7 show that e D 
drops exponentially with n q . This result is similar to 
the one found in [22] for noisy gate errors and can be 
explained by means of the following argument. Static 
imperfections can couple states very far in momentum 
space via a single spin flip. As a consequence, they cre- 
ate n, 



u q peaks [22] with probability W p 



e 2 ff i in each 



peak. Here e e ff <~ Snjjg = en q is the effective perturba- 
tion strength, with r? q T q time between Hadamard gates 
acting on a given qubit. These gates transfer the accumu- 
lated phase error e c s into amplitude errors. Integrating 
the contribution of each peak, one gets 



with 



< (An) 2 >r 



D f 



W P N 2 



e 2 n q N 2 



D e t, 



(13) 



(14) 



One can estimate the critical value eo to double the ex- 
act (e = 0) diffusion coefficient from D t = D n (e = 0), 
giving 



n 2 q N ' 



(15) 




FIG. 5. Dependence of the diffusion coefficient D n on the 
imperfection strength e for k = 2, J = (full symbols), n q — 6 
(triangles down), 8 (stars), 10 (triangles up), 12 (squares), 15 
(diamonds), 18 (circles), and for J = S, n q — 10 (empty tri- 
angles). The straight lines show the theoretical dependence 
D oc e 2 (full line, see Eq. (14)) and the result without imper- 
fections, D(e — 0) « 16 (chain line). 
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FIG. 6. Dependence of the saturation value £oo of the in- 
verse participation ratio on the imperfection strength e, with 
same parameter values and same meaning of symbols as in 
the previous figure. The straight line shows the result with- 
out imperfections, £00 (e = 0) « 14. 

In Fig. 6 we show the dependence of the inverse par- 
ticipation ratio £(e) on e for different number of qubits 
n q . From each curve we extract two critical imperfection 
strengths: 

(i) E^E, to get an inverse participation ratio equal to 
a given fraction of the full Hilbert space, for example 
£ = N/A. This threshold is significant of the transition 
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to ergodic completely delocalized wave functions, 
(ii) ej, to double the exact inverse participation ratio. We 
stress that this quantity gives a rough estimate of the im- 
perfection threshold for reliable quantum computation of 
localization in the absence of error correction. 

The dependence of c^e and on n q are shown in Fig. 
7. These quantities drop polynomially with n q , in sharp 
contrast with the exponential drop of en- This algebraic 
threshold can be understood as follows. The eigenstates 
of the unperturbed (e = 0, J=0) Floquet operator U in 
(2) can be written as 



N 



~(™) 7 



(16) 



where the quantum register (momentum) states. 

In the localized regime, c^'s are randomly fluctuating 
inside the localization domain of size £, and exponentially 
small outside it. Wave function normalization imposes 
\c^\ <~ \j\R-- Due to exponential localization, static 
imperfections couple significantly the unperturbed eigen- 
functions only when their localization domains overlap. 
We estimate in this case the transition matrix elements 
according to perturbation theory. For J = 0, they have 
a typical value 



K yP ~K4 0) lE^^I^ 0) )l 

i=l 

m,n— 1 i—1 

™ q /2 \ ci m) 4 m) V m) | ~ en^r 1 / 2 . 



(17) 



m.—l 



(m) 



In this expression, the typical phase error is b^Jn^r] 
(sum of n q random detunings Si's), with rp m > random 
sign, and T g n g ~ r g n 2 is the time used by the quantum 
computer to simulate one map step. The last estimate in 
(17) results from the sum of order £ terms of amplitude 
|ci"^c^*| <~ l/£ and random phases. Since the spacing 
between significantly coupled quasi-energy eigenstates is 
AE ~ l/£ , the threshold for the breaking of perturbation 
theory can be estimated as 



V t y p /AE 

The analytical result 



e^ E n 5 J 2 V£' 



1. 



5/2 



VI 



(18) 



(19) 



is confirmed by the numerical data of Fig. 7. For the 
case J = 5, the threshold e^E is reduced (see Fig. 4) 
since residual inter-qubit interactions introduce further 



couplings between Floquet eigenstates. However, an es- 
timate similar to (17), which does not modify the func- 
tional dependence (18), can be derived. We stress the 
striking difference between this polynomial scaling and 
the exponential scaling for the mixing of unperturbed 
eigenstates obtained in the ergodic regime (in which 
£ ~ N) and in the more general quasi-intcgrablc regime 
[15,16]. We also note that the different sensitivity of local 
and non local quantities was pointed out in [22]. How- 
ever, the authors of Ref. [22] considered the effect of noisy 
gates, while we consider internal static imperfections. 




FIG. 7. Dependence of the critical imperfection strengths 
on the number of qubits for k — 2, J — 0: thresh- 
olds en (circles), e^E (full diamonds), and (empty di- 
amonds). The full line gives the theoretical dependence 
e^E = Ariq ^ 2 , with the fitting constant A w 0.5. The dashed 
lines gives e D = B^J D{e — Q)2~ nq n q 2 , with the fitting con- 
stant B w 3.6. 



VI. SPECTRAL STATISTICS 

In this Section, we show that spectral statistics is an 
ideally suited tool to detect the destruction of localization 
by static imperfections. We study the spectral statistics 
of the Floquet operator for a quantum computer running 
the quantum sawtooth map algorithm in the presence of 
static imperfections, 



exp 



(ItH(t] 



(20) 



where H(t) is the Hamiltonian (9) and n g the number 
of gates per map iteration. We construct numerically 
the Floquet operator in the computational (momentum) 
basis, using the fact that a single map iteration of each 
quantum register state gives a column in the matrix rep- 
resentation of this operator. Then we diagonalize the 
Floquet matrix and get the so-called quasienergy eigen- 
values Xa^ and eigenvectors (j>a\ 
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CU£>=exp(iA£>)0( 



(21) 



VII. CONCLUSIONS 



A convenient way to characterize the spectral proper- 
ties of the system is to study the level spacing statistics 
P(s), where P(s)ds gives the probability to find two ad- 
jacent levels (quasienergies) whose energy difference, nor- 
malized to the average level spacing, belongs to the in- 
terval [s, s + ds] (see e.g. [23,24]). In the localized regime, 
Floquct eigenvectors with very close eigenvalues may lay 
so far apart that their overlap is negligible. As a conse- 
quence, eigenvalues are uncorrelated, that is their spec- 
tral statistic is given by the Poisson distribution, 



P P {s) = exp(-s). 



(22) 



On the contrary, in the delocalized regime wave func- 
tions are ergodic, and their overlap gives a significant 
coupling matrix element between states nearby in en- 
ergy. In this case the spectral statistics P(s) follows the 
Wigner-Dyson distribution, 



Pwd{s) = — 5-exp 



As' 

7T 



(23) 



typical of random matrices in the absence of time re- 
versal symmetry [23,24] (static imperfections break this 
symmetry). In Fig. 8 we show that static imperfections 
indeed induce a crossover from the localized regime with 
the Poisson statistics to quantum chaos characterized 
by Wigner-Dyson statistics. We have also studied this 
crossover as a function of the number of qubits (data not 
shown): the threshold e c (n q ) for the emergence of quan- 
tum chaos is consistent with the scaling e c (n 9 ) oc n q 5//2 , 
in agreement with the threshold (19) obtained for the 
mixing of unperturbed eigenfunctions. 
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FIG. 8. Level spacing statistics for k — 2, 



. = 11, J = 0, 

e = 10~ 5 (circles) and e = 2.6 x 10~ 3 (squares). The dashed 
and full curves give the Poisson (22) and Wigner-Dyson dis- 
tribution (23), respectively. In order to reduce statistical fluc- 
tuations, data are averaged over No = 5 random realizations 
of Si's, so that the total number of spacings is NdN « 10 4 . 



In summary, we have shown that a quantum computer 
operating with a small number of qubits can simulate ef- 
ficiently quantum localization effects. The evaluation of 
the localization length I with accuracy v requires a num- 
ber of computer runs of order Xjv 2 , followed by a pro- 
jective measurement in the computational (momentum) 
basis. We stress that, in the presence of static imperfec- 
tions, a reliable computation of the localization length is 
possible even without quantum error correction, up to an 
imperfection strength threshold which drops only alge- 
braically with the number of qubits. We also stress that 
localization is a purely quantum phenomenon, which is 
quite fragile in the presence of noise [25,22]. Therefore we 
believe that the simulation of the physics of localization 
can be an interesting testing ground for the coming gen- 
eration of quantum processors operating in the presence 
of decoherence and static imperfections. 
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